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ABSTRACT 


Significant ion and electron flux enhancements immediately 
upstream of the Venus bowshock have been observed by the Electron 
Temperature Probe on the Pioneer Venus Orbiter. It is shown that 
mass loading of the solar wind by oxygen ions accounts for only 
about 10% of the observed effect. 


INTRODUCTION 


Recent observations by the Electron Temperature Probe on the 
Pioneer Venus Orbiter indicate enhanced electron and ion currents 
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immediately upstream of the Venus bowshock (Brace et al., 1985; 
Brace, 1987). Non-turbulent mass-loading of the upstream particle 
flux will be examined in this paper as a potential mechanism 
which could give rise to such a phenomenon. 

Venus is known to have a hot oxygen corona (Nagy et al., 1981; 
Paxton, 1983) resulting in an extended atmosphere. Oxygen is the 
dominant neutral gas species on the dayside out to distances of 
about 3000 km during solar maximum conditions (Kliore et al., 
1986). The bowshock is located at an altitude of only about 2300 
km at the subsolar point during solar maximum (Luhmann, 1986) 
which is inside the region of this hot oxygen corona. The oxygen 
atoms are ionized by either photoionization or charge exchange 
with solar wind protons. The 0+ ions are then picked up by the 
solar wind, thus mass loading the flow. As a result of this energy 
transfer from the solar wind protons to the 0+ ions, the bulk flow 
velocity decreases. This decrease in the solar wind flow velocity 
as well as the addition of 0 + result in an increase in the ion (and 
electron) number density of the mass loaded solar wind. 


DISCUSSION OF THEORY AND UNDERLYING ASSUMPTIONS 

The calculations presented here make use of earlier work on the 
solar wind interaction with comets (Wallis and Ong, 1975; Galeyev 
et al., 1985) with appropriate modifications for the situation at 
Venus. The purpose of this paper is to establish qualitatively the 
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feasibility of mass loading as a cause of the precursor; therefore 
the calculations are simplified as much as possible. The 
calculations described in this paper are based on the following 
assumptions: I) One-dimensional flow, 2) magnetic field 

perpendicular to the solar wind velocity vector, 3) no diffusion in 
either pitch angle or energy of the implanted 0+ ions, and 4) 
conservation of the magnetic moment (first invariant) of the 
implanted ions. 

The photoionization rate of oxygen atoms is given by the 
expression 


Rph * no|ap(v)0(v)dv (1) 

v* 

where ' 

Rph - number of 0+ ions produced by photoionization 

per unit volume per unit time, 
no - hot atomic oxygen density, 

Ej - oxygen ionization potential, 

h - Planck's constant 

o p (v) * photoionization cross section of oxygen for 
photons of frequency v, 

<X>(v)dv - photon flux in the frequency range dv about v. 


Let 


bo 

1 Gp(v)<I>(v)dv, 


( 2 ) 
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then Rph - n 0 I. (3) 

The production rate of 0+ ions by charge exchange is given by 

Rce ■ n o ^cel^pJ^SW (4) 


where 

Rce ■ number of 0 + ions produced by charge exchange 
between O and H + per unit volume per unit time, 
<*ce(Ep) * charge exchange cross section of O with solar 
wind H+ of energy Ep (2 x 10*9 erg for a 
typical solar wind velocity of 450 km/s at 
solar maximum), 

4>SW ■ H+ flux (solar wind) - nooUoo, 
and n^, Uoo are the unperturbed solar wind density and velocity 

respectively. 

Hence, the total 0 + production rate is given by 

Rj * Rph + Rce * not 2 + oce(^p) Oswl- t5) 

Rt is considered to be only a function of the altitude. The two 

cross sections and fluxes are considered to be constant in space 
and time. The altitude dependence of Rj is a result of the altitude 
dependence of n 0 . The calculations presented here are based on 
separate continuity equations, one for the implanted 0 + ions and 
one for the contaminated solar wind, and a momentum equation for 
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the contaminated solar wind. Under the given assumptions the 
upstream solar wind is described by the following equations 

(d/dx)[pjuf(u,p)] = mjn 0 (I + ace^SW) 5 ^ • nriju 2 /2B) (6) 

(d/dx)(pu) - mjn 0 (I + o ce C>sw) ( 7 ) 


(d/dx)(pu 2 + pi + B 2 /8rc) - 0. (8) 

Here the x-axis is along the Venus-sun line with origin at the sun, 
pj is the mass density of the implanted 0 + ions, u is the solar wind 
velocity, f(u,p) is the probability function for magnetic moment of 
the implanted 0 + ions, p-mjvj^/2B is the magnetic moment of the 

implanted 0+ ions, B is the magnitude of the interplanetary 
magnetic field, 5(x) is the Dirac delta function, p is the mass 
density of the mass loaded solar wind, and mj is the mass of the 0+ 
ion. The solution to these equations is (Biermann et al., 1967; 
Galeyev et al., 1985) 

£u - 1 + (mj/pooUoo)(£ + cJce^SW)J no(x')dx' (9) 

p ■ 2pu[1 - (1 - 3pu/4)1/2] t (10) 

u - (2/3pu)[1 + (1 - 3pu/4)1 /2 ] (11) 


Piuf(u.p) 


H-(f 3 -os) 

a iT 

1 r 3 


pooUoo0(Poo - P)0(P - PooU 3 ) 


( 12 ) 



6 


where p, u, and p are the solar wind parameters normalized to their 
unperturbed values p*,, u*,, p«> and where p*, = mjUoo 2 /2Boo and @(x) = 

1 for x > 0 and @(x) = 0 for x < 0. Hence, under the assumption of 
the conservation of the first invariant, the magnetic moment of the 
implanted ions is restricted to the range 

Aq 

PooU^ < P < Poo. 

The total mass density of the implanted ions is therefore 

/* do 

[ pjf(u,p)dp - (1/3)pJj-1(-3 + 4u-1 - u-2). (13) 

RESULTS AND DISCUSSION 

The hot atomic oxygen density profile for solar maximum 
conditions, shown in Figure 1, was obtained from ultraviolet 
spectrometer measurements made onboard the Pioneer Venus 
Orbiter (Paxton, 1983). This results in a column density, defined 
by the integral in equation (9), of 2.4x1011 cm'2. The atomic 
oxygen photoionization rate £ of 4.5x1 0‘^s‘1 at 1 AU for solar 
maximum conditions was obtained from Table 7.19 of Banks and 
Kockarts (1973), which gives a value of 8.6x10"^ s‘1 at the Venus 
orbit of 0.723 AU. The charge exchange cross section for neutral 
oxygen atoms with solar wind protons has been taken to be 
1.5x10-15 cm 2 (equation 9.88 of Banks and Kockarts, 1973). The 
unperturbed solar wind parameters at 1 AU are taken to be n*, = 6 
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cm' 3 , Uoo = 4.5x1 0 7 cm/s which translates into a flux of 5.2x1 0 8 
cm -2 s -1 at 0.723 AU. 


Upon substituting these values into equation (9) one obtains 


A A 

P u 


+ 1.4x10-2 


or 


p u - (1 + 1.4x10-2) pooUoo 


which represents a 1.4% increase in the mass flow. If this solution 
is substituted into equation (10), one obtains 


p - 1 + 3.5x10-2 


which represents a 3.5% increase in the mass density. This 
increase contains the combined effects of the addition of 
implanted ions and the increase in number density of the solar 
wind protons caused by the decrease in the flow velocity. The 
latter, of course, is a consequence of mass loading and can be 
obtained from equation (11), 

u = 1 - 2.1 x 10"2. 


Hence, the mass loading results in a slowdown of 2.1% in the solar 
wind speed in front of the Venus bowshock. We are now in a 
position to calculate the mass density of the implanted 0+ ions and 
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the increase in the number density of H+ ions. Using the 
parameters listed above, the mass density of implanted 0+ ions is 
obtained from equation (13) 

J-io 

pj\f(u,p)dp - 1.4 x10" 2 poo. 

USfd © 

Thus a 1 .4% mass density increase is due to the addition of 0 + to 
the flow. The remaining 2.1% increase (to give a total mass 
density increase of 3.5%) is due to the density increase in solar 
wind protons. 

These results are of course based on the simplifying assumptions 
stated in the preceding section. Recently a number of papers was 
published presenting more general theories (e.g. Galeyev, 1986; 
Gombosi, 1988) than the one on which this paper is based. In 
particular, these authors include pitch angle scattering and thus 
the violation of the first invariant. In that case the pickup ions are 
distributed on a pickup shell in velocity space instead of on a 
pickup ring as they are in this treatment. Ion pickup is obviously 
maximized when the ions travel at right angles to the magnetic 
field. 

A second effect of pitch angle scattering is a change in the value 
of the specific heat ratio y which has been taken to be 2 in the 
calculations presented here. If pitch angle scattering is included, y 
should be taken as 5/3 (appropriate for three degrees of freedom). 
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The effect of a decreased value of y, can be obtained from the 
equation for the deceleration of the mass-loaded solar wind 

du/dx = (y+ 1)u 2 pn/2(yp - pu 2 )* 

(Wallis and Ong, 1985). It is easily seen that du/dx decreases with 
increasing y or - in other words - the solar wind deceleration 
increases with increasing y. Hence, the reduction in the value of y 
leads to a reduced solar wind deceleration which in turn implies a 
reduced plasma density increase. The inclusion of pitch angle 
scattering in our calculations would have resulted in a smaller 
density enhancement for two reasons, a redistribution of the 
pickup ions to pitch angles other than 90° and a reduction in the 
specific heat ratio y. Our results therefore represent an upper 
bound of the electron density enhancement due to mass loading of 
the upstream solar wind by 0+ ions. 

The description of the pickup ions (equation 6) is based on the 
assumption of a small gyroradius compared to the dimensions of 
the interaction region under consideration. In the present case this 
is the thickness of the precursor. The gyroradius of 0+ ions of near 
solar wind speed perpendicular to the IMF is of the order of several 
thousand km which is about the magnitude of the precursor 
thickness. Phillips et al. (1987) have recently shown that finite 
gyroradius effects are significant in the Venus magnetosheath due 
to the cylindrical asymmetry of the motional electric field in the 
magnetosheath between the two Venus hemispheres. In the 
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unshocked solar wind, finite gyroradius effects will result in a 
redistribution of the pickup ions. The mass leading locations, on 
the other hand, are distributed according to the locations where 
the oxygen ionizations occur which in turn follow the distribution 
of neutral oxygen at a given radial distance from Venus. The solar 
wind decelerating force field acts at the mass loading locations 
and therefore also follows the neutral oxygen distribution. 

Therefore the plasma density increase due to mass loading is not 
affected by finite gyroradius effects, and, as shown above, this is 
the major part of the total density increase. 

CONCLUSIONS 

The calculations presented here show that non-turbulent mass 
loading of the solar wind by 0 + will result in a number density 
increase of at most 2.2%. If the observed 30% electron flux 
enhancement in the precursor is assumed to be entirely due to a 
corresponding density enhancement in the upstream solar wind, it 
is of interest to calculate the minimum required atomic oxygen 
column density which is necessary to obtain such a density change. 
This turns out to be a column density of 3.2x1 0^2 cm"2 which is 
about one order of magnitude larger than the column density used 
in the above calculations. No atomic oxygen column density of that 
magnitude has been observed by Pioneer Venus (Paxton, personal 
communication). Therefore these calculations indicate that non- 
turbulent mass loading of the solar wind flow in front of the 
bowshock cannot be the sole mechanism responsible for the 


observed enhancements in the electron and ion fluxes. Sagdeyev et 
al. (1986) have recently shown that mass loading of the solar wind 
to distances beyond the cometary bowshock excites Alfven 
turbulence the amplitude of which is of the order of the 
interplanetary magnetic field. Such turbulence has indeed been 
observed near comet Halley (including the region outside its 
bowshock) not only in the magnetic field but also in the electron 
density data (Riedler et al., 1986; Neubauer et al., 1986; Neubauer, 
1986; Reme et al.. 1986a, b). These observations also include an 
electron density enhancement upstream of the cometary bowshock. 
Therefore the observed increases in the Langmuir probe collection 
current upstream of the Venus bowshock must either be the result 
of some other physical process in the preshocked solar wind flow, 
like the above mentioned Alfven turbulence, or the result of some 
unexpected instrumental and/or spacecraft phenomenon. 
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FIGURE CAPTION 
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FIGURE 1 


: Hot atomic oxygen density distribution in the upper 
atmosphere of Venus (Figure B-8 of Paxton, 1983) 
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